Non-Markovian qubit dynamics in the presence of 1/f noise 
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presence of 1/f noise, without Markov approximation. We show that the non-Markovian qubit time- 
evolution exhibits asymmetries and beatings that can be observed experimentally and cannot be 
explained within a Markovian theory. The present theory for 1/f noise is relevant for both spin- and 
superconducting qubit realizations in solid-state devices, where 1/f noise is ubiquitous. 
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I. INTRODUCTION 

Random telegraph noise has been encountered in a 
wide range of situations in many different areas of 
physicsA. A typical example in condensed matter physics 
is that of a resistor coupled to an ensemble of randomly 
switching impurities, producing voltage fluctuations with 
a spectral density that scales inversely proportional with 
the frequency, hence the name "1/f noise". The quest to 
build and coherently control quantum two-level systems 
functioning as qubits in various solid state systems has 
once more highlighted the importance of understanding 
1/f noise, being a limitation to the quantum coherence 
of such devices. 

The description of low- frequency noise (such as 1/f 
noise) is complicated by the presence of long-time corre- 
lations in the fluctuating environment which prohibit the 
use of the Markov approximation. Only in few cases, non- 
Markovian effects have been taken into account exactly, 
e.g., for the relaxation of an atom to thermal cqulibrium 2 . 
Here, we are interested in the decoherence and relaxation 
of a qubit, i.e., a single two-level system (spin 1/2). For 
the spin-boson model, i.e., a qubit coupled to a bath of 
harmonic oscillators, the dynamics has been calculated 
within a rigorous Born approximation without making a 
Markov approximation 3 '^. Here, we carry out a similar 
analysis for 1/f noise and find even stronger effects than 
in the spin-boson case (see Fig.Q]). 

Charge and to some extent (via the spin-orbit 
interaction) spin qubits in quantum dots^ formed 
in semiconductor 6 or carbon 7 structures are sub- 
ject to 1/f noise. In superconducting (SC) Joseph- 
son junctions, SC interference devices (SQUIDs), and 
SC qubits, 1/f noise has been extensively studied 
experimentally^ii&iiii^iii^ and theoretically^ 6 .^. 

Even where the origin of 1/f noise is known, the in- 
duced decoherence is not fully understood. Most theoret- 
ical work is either restricted to longitudinal fluctuations 
or employs a Markov approximation. Here, we present a 
calculation of the qubit dynamics in the presence of 1/f 
noise which is exact within the lowest-order Born approx- 
imation. In particular, we make no use of a Markov ap- 
proximation. In contrast to earlier calculation a 18 ' 19 i 21 i 22 , 
we allow for arbitrary qubit Hamiltonians and include 
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FIG. 1: (Color online) Non-Markovian time-evolution of 
the unbiased (e = 0) qubit (spin) z-component {o~ z (t)), for 
A/A 2 = 0.05 and 70/A = 0.05 (solid black line). The 
Markovian pole contribution (<Xz(i)) P oies is plotted as a dashed 
line for comparison. The essential non-Markovian part is 
non-exponential and given by the branch cut contribution 
{o- z {t)) hc (red solid line). Inset: Plot for A/ A' 2 = 0.005 and 
70/A = 10~ 10 . Here, the essential non-Markovian part is the 
long-time asymmetry which carries information about the ini- 
tial state. 



transverse as well as longitudinal (phase) 1/f noise. Non- 
Gaussian 1/f noise originating from few fluctuators was 
studied i n 23 i 24 ' 25 , while numerical studies using an adia- 
batic approximation were carried out in^ 6 -. The coupling 
to a single fluctuator was also studied^ 7 .. 



II. MODEL 

We model the qubit (spin 1/2) coupled to a bath of 
two-level fluctuators with the Hamiltonian 



H = H S + H E 



SB 



with 



Hsb 



(1) 

(2) 
(3) 



where a x and a z are Pauli matrices describing the qubit 
and X = Yli=i v i a l where a\ operates on the z-th fluc- 
tuator. In a SC qubit, A and e denote the tunneling 



2 



and energy bias between the two qubit states. In a 
spin qubit, e is the Zeeman splitting and A a transverse 
field. The bath Hamiltonian Hb need not be provided 
explicitly; it is sufficient to know the auto-correlator 
C(t) = (X(0)X(t)) of the bath operator X(t), where 
(. . .) = Ttb(---Ab) denotes a trace over the bath de- 
grees of freedom with the bath density matrix ps- We 
can further assume that the fluctuators are unbiased, 
(X(t)) — 0. For independent two- level fluctuators with 
switching rates ji, one obtains 

c{t) =x>ifo(*) ff '(°)> = E w * e " 74|t| - w 

i i 

The noise spectral density is the Fourier transform 

/oo 
dtC(t)e- iut = J2(^H)/( 7 f + to 2 ). (5) 
-oo i 

While this correlator describes essentially classical bath 
dynamics (as is commonly assumed for 1/f noise), it 
should be emphasized that our model is not classical, 
because the [Hsb, Hs] ^ 0. In the case of a large num- 
ber of fluctuators, the sum in C{t) can be converted into 
an integral. For 1/f noise, one typically assumes a distri- 
bution of fluctuators of the form P(v, 7) cx l/'yv", where 
both v and 7 are limited by upper and lower cut-offs 31 . 
The spectral density of the ensemble of fluctuators then 
becomes 



S(u>) tx 



7c 2v 2 j 
dvdjP(v,-f) — 

70 



r 



(6) 



For 70 = this yields 1/f noise of the form S(u>) cx l/|w|. 
The divergence at low frequencies is cut off by the finite 
duration of a qubit measurement, if not by other effects 
at even shorter times. A low-frequency cut-off 70 > 
yields 



S(uj) = 2txA 



arctan(u;/7o) 1 



(7) 



where A depends on the cut-offs and the exponent (3. For 
70 — > 0, we recover S{uj) — > 2~kA/\uj\. Inverting the above 
Fourier transform, we obtain 

C(t) = -AEi(- 7o \t\), (8) 

where Ei denotes the exponential integral function. 

III. QUBIT DYNAMICS 



where the self-energy superoperator £(t) gives rise to 
memory effects, i.e., the time evolution of ps{t) depends 
on the state ps(t') at all earlier times t' < t. There- 
fore, the qubit dynamics is inherently non-Markovian. 
Expanding the right-hand side of the GME in orders 
of Hsb and only keeping the lowest (second) order, 
one obtains £ in (lowest-order) Born approximation 
E(t)p s - -iTr B [H SB ,e- itH "[H SB , p s ® Ps]e itff °], where 
H = Hs + Hb ■ 

Introducing the Bloch vector (<r(t)) = Trserps(i), 
where <x = (cr x , a y , a z ) is a vector of Pauli operators, 
we write the GME as a generalized Bloch equation 

(&) = R * (a) + k, (10) 

where the star denotes convolution and 3,4 



-£ r i(*) 



sS(t) + %K+(t) 







R(t) = I e5(t) - %K+(t) -T y (t) -AS(t) 







AS(t) 







(11) 



with E = VA 2 + e 2 and 3 - 4 T^t) 
(2A/E) 2 cos{Et)C'(t), T y (t) = (2A/E) 2 {1 + 
(e/A) 2 cos(Et))C'(t), and K+{t) 

(4eA/£: 2 )sin(M)C"(i), where C'(t) and C"{t) de- 
note the real and imaginary parts of C(t). Since for 1/f 
noise, C"(t) = 0, we find k(t) = 0M. As shown in 3 .^, 
Eq. (fTU)) can be solved by means of the Laplace transform 
(LT) f(s) = f™f(t)e- ts dt, where 

(<T( S ))=( S -R( S ))- 1 ({cT(t = 0))-k( S )). (12) 

The LT R(s) of R(t), has entries according to Eq. (JTTJ, 
with 6(t) replaced by 1, and, for 1/f noise 

r x (s) = {2A/E 2 )A 2 (C{s + iE) + C(s-iE)), (13) 
T y (s) = {2A/E 2 )(2A 2 C{s) 

+e 2 (C(s + iE) + C(s-iE))) , (14) 
K+(s) = i(2A/E 2 )Ae(C{s + iE)-C(s-iE)), (15) 

where the LT of the correlator C(t) in Eq. ([4]) is 



C( S ) = -log(l + S /7o) 



(16) 



We recover (<r(t)) from (<r(s)) by way of an inverse LT as 
carried out below, first for the special case of an unbiased 
qubit (e = 0) and then for the general case. 



The density matrix p of the total system, consisting 
of the qubit and the bath, obeys the Liouville equation, 
p(t) = — i[H, pit)]. The time evolution of the reduced 
density matrix of the qubit alone ps(t) — TtbP is then 
determined by the generalized master equation (GME) 3,4 

p s (t) = -i[H s , ps(t)] -if E(t - t')p s (t')dt\ (9) 

Jo 



IV. UNBIASED QUBIT 

We first assume that the qubit is prepared at time 
t = in one of the eigenstates |0) = | f) of <r z , i.e., 
(er) = (0,0,1), and that the qubit is unbiased, e = 0. 
If the fluctuators were absent the qubit would undergo 
a precession about the x axis, (c z (t)) = cos(Af). Due 
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to the presence of the fluctuators, we find (see also Ap- 
pendix [A} 



s 2 +4Alog(l + s/ 7 o) 
s(s 2 + A 2 +4^1og(l + s/ 7o ))' 



(17) 



We expand (a z (s)) in leading order of A, 



<*»oo> 



A 2 



iAA2 Ww +0{A2y (18) 



The coherent spin oscillations in the time do- 
main are obtained from the inverse LT, the so- 
called Bromwich integrate (see Fig. [2]), (a z (t)) = 

lim^ i0 f^Z l TI Wz(s)}e t8 ds - The integral contour can 
be closed in the left complex half-plane Re(s) < 
(Fig. [J). The behavior of (<r z (t)) is therefore given by 
the analytic structure of (<J z (s)) in the left half-plane, 
see Fig. [2l In the absence of the fluctuating environment 
(A = 0), (<7 z (s)) has two poles at s = ±iA which yield 
(&z(t)) = cos(Af), as expected. The coupling to the en- 
vironment has two effects: (i) a shift of the poles, and 
(ii) the appearance of a branch point (bp) due to the 
logarithm in Eq. (|18p and the associated branch cut (be) 
that we choose to lie on the real axis between —70 and 
—00. Here, it should be noted that in the case of an 
unbiased qubit, the presence of 1/f noise does not lead 
to the appearance of a pole on the real axis, and thus 
there is only pure dephasing and no T\ type decay (spin 
relaxation), in contrast to other types of environment 4 . 
The exact shift of the poles has been calculated numer- 
ically from Eq. (JTTJ). To lowest order in A, we find 



A r = A' r + iA'; 



A 



1 



A- 
To 



±2i£ arctan— , 

A 70 ' 



where the real part A' r is the renormalizcd frequency of 
the coherent oscillations, while the imaginary part A" 
describes an exponential decay of those oscillations. If a 
Markovian approximation were made by setting s = in 
Ti(s), T y (s), and K+(s), then the be would be missed 
completely and only an exponential decay with a rate 
2 A/70 would be obtained. The Markov approximation 
is only justified if 70 S> A, i.e., if the bath dynamics is 
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FIG. 2: (Color online) Analytic structure of (a z (s)) in the 
complex s plane, for (a) the unbiased case, e — and (b) the 
biased case, e / 0. Red dots denote poles, blue lines branch 
cuts. 



much faster than the system dynamics. Here, we entirely 
avoid making a Markov approximation. 

The Bromwich integral can then be divided into 
two parts, (a z (t)) = (cr z (t)) poles + (a z (t)) hc . The 
integration in the first term along the contour C, 
not including the line integrals along the be (Fig. ^ 
yields the sums of the residues from the poles 
MOW = ^-J c ds{a z {s))e st = r' cos(A;t) e - A " t - 
r"sin(A;t)e- A "', where r' = 1 - (2A/A 2 ) log(l + 
A 2 /7o) + 0(A 2 ) and r" = (4A/A 2 ) arctan(A/7 ) + 
0(A 2 ). For A = 0, this reduces to cos(At). 

The branch-cut contribution to lowest order in A is 



AA 

(°-,(t))bc = ^i(7o/A,At) 



(19) 



with the integral I n {a,b) = J a °° dy jy „ ( e y2 + 1) 2 , where we 
have used Eq. (|18p and introduced dimensionless vari- 
ables and where a > and b > 0. For n = 1, we find 
(Fig.© 



h(a,b) = ^Rc[{ib + 2)e~ tb {-iK + Vi{ib- ab))} 



1 1 



2 1 + (X s 



e - ab -Ei(-aft). 



(20) 



For a = 70/A > 1 and b > (t > 0), the effect of 
the environment from the be integral is exponentially 
suppressed: Ii(a,b) < e~ ah ' jb and thus |(cr z (f))bc| < 
{AA/ A 3 i)e~ 70 * . The physically more interesting regime 
is a = 70/A <C 1. Within this regime, we can distinguish 
two temporal regimes: short times ab <C 1 (t <C "Jq 1 ) and 
long times ab ^> 1 (t ^> 7 ~ 1 )- In the short-time case, the 
integral is cut off from above by a combination of the 
y~ 5 and the exponential factor. The effect of the latter 
can be approximated by cutting off the integral at 1/6, 
with the result I\{a,b) w — 7i(l/&,0) + /i(a,0), where 
h(a,0) = -i(l + a 2 )" 1 + ilog(l + a -2 ) is the be integral 
for t = (b = 0). Note that h(a,0) > due to the log- 
arithmic term. In the long-time case, the integral is cut 
off by the exponential whereas the (y 2 + l) 2 factor in the 
denominator becomes irrelevant, Ii(a,b) w — Ei(— ab). 

At this point, the parameter that controls the strength 
of the non-Markovian effects due to 1 /f noise can be iden- 
tified as f = (A/A 2 ) log(l + A 2 /^l). The regime of valid- 
ity of the Born approximation (the only approximation 
required in this paper) is confined by the condition £ <C 1. 
The resulting damped qubit oscillation is plotted in Fig.[T] 
for A/ A 2 = 0.05 and 70/A = 0.05 where £ « 0.1. If the 
infrared cutoff is lowered, the non-Markovian effects due 
to 1/f noise become more pronounced. However, since 
the dependence on the infrared cutoff 70 is only logarith- 
mic, the result does not change drastically even if 70 is 
much smaller than in our example, as long as A is cho- 
sen sufficiently small to ensure the validity of the Born 
approximation. E.g., for A w 10 GHz and 70 « 1 Hz (cf. 
Ref. [l3) then 7o /A = lO" 10 . With^ A/ A 2 = 0.005, 
one finds a long-lived asymmetry as shown in the inset of 
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Fig. ([T]) . The intermediate asymptotics of this contribu- 
tion is (cr z )bc ~ £ ~ 0-1 j while for longer times this con- 
tribution also decays logarithmically to zero. A similar 
long-time behavior has been found also for longitudinal 
coupling 10 . 



THE BIASED CASE 



We again assume that the qubit prepared at time 
t = in one of the eigenstates |0) = | T) of a z , i.e., 
(tr) = (0,0,1), but now the qubit is biased, e ^ 0. In 
the absence of the fluctuators (A = 0), the qubit would 
now undergo a precession about an axis in the xz plane 
with frequency E/2ir, where E = \/ A 2 + e 2 . In this un- 
perturbed situation, (<j z (s)) has three poles at s = ±iE 
and s — 0, the former two giving rise to undamped os- 
cillations of (a z (t)) with frequency E/2ir and amplitude 
A 2 /E 2 , while the latter allows for a non- vanishing sta- 
tionary value e 2 /E 2 of (<r z (t)) in the long-time limit. 

Including 1/f noise we find in leading order in A (see 
Appendix |A"|) . 



2 i 2 
S + 6 



A 2 
4A—rRc 



s{s 2 + E 2 ) E 2 



Ms) 



s 2 {s + iE) 



;C{ S + iE) 



{E 2 + s 2 ) 
+ 0(A 2 ). (21) 



Analogously to the unbiased case, the poles are shifted in 
the presence of the fluctuators. In leading order in A, we 
find three poles at = - (4AA 2 / E 3 ) &rcta,n(E / j ) , 

and ±iE r = ±iE ± {iAA 2 /E 3 ) log(l + E 2 ^ 2 ) - 
(2AA 2 /E 3 ) arctan(i?/7o). From the shift of these poles 
(Fig. [2b), we obtain (cr z (t)) polcs = j£ cos( J B^)e- E "* + 

j^e~ 2Er *. However, while in the unbiased case a Marko- 
vian treatment at least qualitatively describes the pole 
contribution correctly, in the biased case, there is an- 
other effect that is elusive in a Markovian analysis. As 
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FIG. 3: (Color online) Branch cut integral function Ii(a,b) 
(solid red line) and two asymptotes (black dashed and blue 
dotted lines). 
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FIG. 4: (Color online) Shift and splitting of the poles of 
(c z (s)) for the biased system (e = 0.3 and 70 = 0.05). Shown 
is the pole located at s = iE for the undamped system 
(A = 0), indicated as a red dot (see Inset b). The pole at 
s — —iE behaves similarly. With increasing A the pole shifts 
toward the vicinity of the branch point (bp), where a second 
pole (orange dot) appears. Shown as red and orange dots are 
the two poles for A = 0.05. The splitting of the poles leads to 
a beating in {tJ z (t)) po \ CB , see Fig. [5] Inset: Analytic structure 
of (ctz(s)) for e / 0, where red dots are poles, and blue lines 
are branch cuts (see Fig. [2}. 



shown in Fig. [2Jd, there are three bp's in the biased case, 
lying at —70 and —70 ± iE. We find that as the two 
poles near ±iE approach the bp's at —70 ± iE as A is 
increased, these poles split into two poles. This behavior 
is illustrated in Fig. [4] The significance of this splitting 
is that it leads to beating patterns already in the pole 
part of (<7 z (t)), as shown in Fig. O It should be noted 
that, again, the precise value of 70 is not critical for the 
possibility to observe the effect, since 70 only enters in 
the argument of a logarithm; even a much smaller value 
of 70 can thus be compensated by only a slight increase 
of the system-environment coupling constant A. 

The three bc's give rise to a contribution to (a z (t)), 



(c z (t))bc = - 



AAA 2 



E* 
■2 



] cos(Et) 



E 2 



E 2 



(sm(Et)I 2 - cos(Et)h) 



(22) 



where the functions /„ are as defined above and are eval- 
uated at the arguments a = ja/E and b = Et. For the 
unbiased case e = and E = A, one retrieves the pre- 
vious result. The integrals I2 and J3 can be calculated 
in closed form, but will not be given here. The damped 
oscillations (a z (t)), consisting of both pole and be con- 
tributions, are plotted in Fig. [5] 
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FIG. 5: (Color online) Oscillation (er z (t)) of the biased qubit 
for e/A = 0.3, A/A 2 = 0.05 and 7o /A = 0.05. The beating 
due to the splitting of the poles at ±iE can be observed in 

(°* W) poles- 



VI. COMPARISON WITH AN EXACTLY 
SOLVABLE CASE 



The circumstance that in the case A = the coupling 
Hamiltonian between the system and the environment 
Hsb commutes with the system Hamiltonian Hs makes 
this special case exactly solvabl o 18 i 19 ' 21 i 22 . A state pre- 
pared transverse to the common direction of the fixed 
precession axis and the fluctuating field, e.g., as (<r(t = 
0)) = (1, 0, 0), for low-frequency noise essentially leads to 
a Gaussian decay behavior (er x (t)) = cos(ei) exp (— ct 2 ). 
The Born approximation which we have employed here 
can only be expected to yield this result in lowest-order 
of the coupling constant, i.e., 



.(*)) 



i(et) (1 - ct 2 + 0(c 2 t 4 )) 



(23) 



Here, we show that our result indeed has this form in the 
special case A = 0. 

To this end, we take the limit A — > in the propagator, 
Eq. (|12p. as shown in the Appendix El We then find 



(<7x(s)) = Pxx(s) = 



s + T v (s) 



( S + r y ( S )r+[e-K+(s) 



(24) 



From Eq. (j!6[) and omitting logarithmic corrections, we 
can use C(s) ~ A/s, and thus T y (s) ~ 4As/(s 2 + e 2 ) and 
K+(s) ~ 4Ae/(s 2 + e 2 ). Substituting this into Eq. 
and expanding to lowest order in A, we find 



(c*0)} 



,s z + e z 



A 



s(3e 2 



s 2 ) 



(s 2 + e 2 ) 3 



(25) 



which equals the LT of Eq. I|23p to lowest order, with the 
identification c = Aj 2. Therefore, our result is consis- 
tent with the known exact result for A = 0, but, within 
the Born approximation, goes far beyond it, in that it 
includes arbitrary values of e and A. 



VII. DISCUSSION 

We find the following essentially non-Markovian fea- 
tures in the decay of the z-component of the spin: (i) The 
spin decay is non-exponential and asymmetric. For rela- 
tively large infrared cutoff 70, there is an "initial loss" of 
coherence on a typical time scale I/70, as seen in Figs.Q] 
and[5l More importantly, for the typical case of small 70, 
there is a long-time asymmetry favouring the qubit near 
its initial state, (ii) In the biased case, 1/f noise can lead 
to a two-frequency oscillation, exhibiting a characteris- 
tic beating pattern. Here, we have concentrated on the 
longitudinal component (a z (t)) of the qubit under the 
influence of both longitudinal and transverse 1/f noise. 
The transverse component (o~ x (t)) shows similar behav- 
ior. The predicted non-Markovian effects are observable 
in free induction decay (Ramsey fringe) experiments. In- 
deed, such asymmetries are clearly visible in supercon- 
ducting qubit o 10 ' 29 Measurements on a superconducting 
flux qubit have shown deviations from the exponential 
decay and beatings^. The question whether these ef- 
fects are due to the mechanisms described here or not 
require further investigation. 

Acknowledgments. Financial support for this work 
from the Swiss SNF (contract PP002-106310/1) and from 
German DFG SPP 1285 "spintronics" and FOR 912 is 
gratefully acknowledged. 



APPENDIX A: FORM OF THE PROPAGATOR 

The propagator (resolvent) for solving the generalized 
Bloch equation in Laplace space is defined in Eq. (fT2"|) as 



P(s) = (s-R(s)Y 



(Al) 



Using the form of the relaxation matrix i?(s), we obtain 
the following expressions for the matrix elements of P(s), 



Pxx (5) 

Pyy( s ) 

Pzz(s) 
Pxy(s) 
Pxz(s) 

Pyz(s) 



1 / A 2 

W)\ s + Tyis) + ~ 

1 / E 2 

---P (s) 
9 r yy\ b )i 



s s 
-Py*(s) 

Pzx{s) = 



(A2) 
(A3) 
(A4) 



A 



Pxy(s), 



~Pzy{ s ) — Pyy{ s )> 



(A6) 
(A7) 



with the definition 



D(s) 



s + r v (s) + 



(A8) 
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The solution in Laplace space is now obtained according 
to Eq. (O, with k = 0, 



(A9) 



j=x,y,z 



E.g., for (<r(t = 0)) = (0,0,1), we find (^(s)) = P iz (s). 
Using Eqs. (|A4]l . and (32J), we recover the known 

results from Ref. [1 in the special case k = 0. The remain- 
ing matrix elements, Eqs. (|A2|) . (|A3|) . and (|A5|) . allow us 
the use different initial conditions. 



1. The case e = 

For an unbiased qubit, e = and thus E — A, so that 
the quantities discussed above are reduced to the form 

D(s) = (s + T y (s) + — )( s + r 1 (s)), (A10) 



(All) 



Pyy(s) = £J c^= (« + r„(«) + ^) ,(A12) 

p 2Z ( s ) = ( s + r y ( s ))p„( s )/ s , (Ai3) 

P y ,(«) - -P zy (s) = -AP ra (s)/s, (A14) 

P* y (s) - Pyx(s) = P xz {s) = P zx (s) = 0. (A15) 



2. The case A = 



In the case of a diagonal system Hamiltonian Hs, we 
set A = and thus E = e, and 



A 2J 



r ^ s ) = ^r 1 ( S ) = 2A(C( S + le ) + C( S -ze)) : 



(A16) 



K+(s) ee -K+(s) = 2iA(C{s + ie)-C(s-ie)), 

(A17) 

D( S ) = (s + T y ( S )) 2 +[e-K+( S )y . (A18) 
With Eqs. (jAlHATf . we obtain 



Pxx{ s ) — Pyy{ s ) ~ 

PM = - 



S + Ty(s) 



(A19) 
(A20) 



e-K+(s) 

P X y{ S ) = -Py X (s) = q^lA, (A21) 

P xz {s) = P zx {s) = P yz (s) = P zy (s) = 0. (A22) 
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